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ABSTRACT 

The problem of small perturbation potential supersonic flow 
around complex configurations is considered. This problem requires 
the solution of an integral equation relating the values of the 
potential on the surface of the body to the values of the normal 
derivative , which is known from the small perturbation boundary 
conditions. The surface of the body is divided into small (hyper- 
boioidal quadrilateral) surface elements, , which are described 
in terms of the Cartesian components of the four corner points. 

The values of the potential (and its normal derivative) within 
each element is assumed to be constant and equal to its value 
at the centroid of the element. This yields a set of linear 
algebraic equations. The coefficients of the equation are given 
by source and doublet integrals over the surface elements, . 

Closed form evaluations of the integrals are presented. 
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SECTION I 

FORMULATION OF THE PROBLEM 


1,1 Introduction 

A general theory for compressible unsteady potential 
aerodynamic flow around lifting bodies having arbitrary shapes 
and motions is given in Refs. 1 and 2. Applications to wings 
in subsonic flows are given in Refs. 3, 4 and 5. A finite 
element formulation for complex configurations in subsonic 
flows is given in Ref. 6. Here the formulation is extended 
to steady supersonic flows. 

The equation of the aerodynamic potential is given by 




* O 


( 1 . 1 ) 


with boundary conditions ( h b is the outwardly directed normal 
to the surface of the body, C b ) 


2-? * - r\ x (on <*fc) 


( 1 . 2 ) 


The Green function for Eq. (1.1) is , -for M > | 


G 


2tt x 


(1.3) 


where H is the Heaviside function and 

e> 'Jm"-- i 


(1.5) 
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The Green theorem for Eq. (1.1) is (Ref. 1) 
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where v,. is the gradient in the Physical Space variables 
(x,y,z), O' is a suitable surface (see Section 2), defined by 
the equation 

S (K,i,e) ~ © 

? (1.7) 

and the function E is given by 

f »■ O Inside <y 

E s i ovtiide <5 

Introducing the supersonic Prandtl-Glauert nondimensional 

variables (■£ is a characteristic length of the body) 


( 1 . 8 ) 
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yields 
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where is the gradient in the Prandtl-Glauert variables# £ 
is the surface of the space (X,Y,Z) defined by the equation 

s cex,r.z) =o 


(1.11) 
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and 


R * J (x-x.)*- (r- r.) -(z-z.f 

Combining Eqs. (1.6) and (1.10) yields 

2tr* A , ,6;(-.22.±±+2l-S£- T ±L j£.) 
'*• “l ax ax ar ?Y az az / 


( 1 . 12 ) 
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Note that 
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(1.13) 


(1.14) 


is the normal to the deformed surface £ . In order to use 
compact vector notations/ it is usual (see, for instance/ Ref. 

7) to introduce the concept of the conormal 

N C = N, i - N t J - A/ t K (1.15) 


and the cogradient 


7* . -2- i 

V T>% 


® 7 "» 

n * ' a* K 


(1.16) 


With this notation, Eq. (1.13) reduces to 


iire +-‘#(T 7 


(1.17) 


where the conormal derivative is given by 

V * ' 7 fc (1.18) 

Equation (1.16) is in agreement with Eq. (6.89) of Ref. 1 , 
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where the opposite convention 01 the direction of the normal 
is used (here the normal n is directed from the region E = 0 
to the region E = 1) . 

On the other hand, the boundary conditions on the surface 
of the body, £* , can be written as 






U»j? I gS> 2 j»24 2li *$. \ I 

Jl Q 0* *X 0X -ay 2Z aZ/ * 3 vxj 


y» (_ 25 *5 _»$_ 25 . 24 . 4 (. I i ) 25 . 24 . 4 . 1 

~ Q l dX vY »Y »Y « 32 V 0* * v »x *X J 

_ 0 (1.19) 



( 1 . 20 ) 


Neglecting terms of the same order of magnitude as those 
neglected in linearizing the equation of the aerodynamic 
potential yields 


Z 4> 


*~2 NJ A 


^ i. ^ on the body (1.21) 

Finally, note that the linearized Bernoulli theorem for 
potential flow yields ■ 


c ' l, s -4 *± 


(1.22) 



or, in Prandtl Glauert variables, 
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(1.23) 


It may be worth noting that the exact boundary conditions 
may be rewritten as 



(1.24) 


1.2 Supersonic Vector Algebra 

The use of the conormal and the cogradient was found to 
be quite cumbersome for the extension of the subsonic finite- 
element formulation to supersonic flow. The algebraic mani- 
pulations become much simplified if a special vector algebra 
is introduced. This algebra is called here, supersonic vector 
algebra or super-algebra . The sum, dot product and cross 
product are defined in the usual way. In addition, it is 
convenient to introduce the supersonic dot product or super- 
product as 

a O i - dn b* - b 1 - b t (1.25) 


With this notation. 


where 





(1.26) 


(1.27) 
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Similarly 

= fj « v (1 - 28 > 

Furthermore, in addition to the usual norm of a vector 

1*1* J1 PIT" (1.29) 

it is convenient to use the supersonic norm or super-norm 

* /l| S, © tf/| (1.30) 

Note that 


5 




(1.31) 


that is , Qo(ji>O^QOQ<Oi Q & m - O ^ if the angle between 
the vector a and the x-axis is less than (greater than; equal 
to) 45°, which implies that the vector a is pointed inside 
(outside; on) the Mach cone (see Fig. 1) . Further development 
of the superalgebra is given in Appendix A, where the first 
super-rule 

(a* l)o(c*B ) « (ao c)(lo3)-(ao3)(l°c) (1<32) 

and the second super-rule 

( 5o a) [(I * c.)o ( 5 1 c)]- (a • l*c) 

& (aoc){(l*c)o(lfa)}+ (ao},)[(ct2)o(ct aj) 

(1.33) 

and the third super-rule 

are derived. ORIGINAL PAGE IS 

OF POOR QUALITY 
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With this notation, Eq. (1.13) may be rewritten as* 




(1.35) 


1.3 Finite Element Formulation 

Assume that the surface X is divided into N small finite 
elements and consider the simplest finite element repre- 
sentation, that is, assume ^ and 

♦* - fi . - ijovf a - 36) 

constant within each element. Equation (1.35) then reduces to 

E.4. - -,l H jj- - i <b. _L((i}.^_±L_<)2- 

*•' T * JToT ' :ir)J JTTT ' 


- i f "" s - - z d> . p 

it, T > in ”* ^ 


where and are the source integral 


(1.37) 


5 * 27T [) 

and the doublet integral 






(1.38) 


(1.39) 


*Note the analogy with the subsonic Green Theorem 
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evaluated on the ith element, 1 ' M . 

In Sections 2 and 3, the evaluation of S and D is derived 
for a hyperboloidal element (see next subsection) . Note that 
for a planar element* 



A/o 7. >5 


(1.40) 


This relation is extremely useful since it is possible to 
evaluate D (which in general would involve the use of the 
finite part of the integral) as the conormal derivative of S 
(which does not contain infinite part). 


1,4 Hyperboloidal Element 
Consider the equations 

X * X. 7,f7 

l f * y. - i’i * i. 7 * 1>17 

i • * J.f ♦ b? U-41) 

or, in vector notations 

p * p< * p,-?- ?,■)* F,V tl - 42) 

This represents a hyperboloid. The lines tj = const and ^ * 
const are clearly straight lines. Consider the hyperboloidal 
element (Fig. 2) defined by the above equation with 

-I £ ? £ I 

(1.43) 

-I £ 9 1 I 

Note that V { (p- P,) s (p'f,) where y is the gradient for 
the variables ( ) 
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The centroid of the element is p i'f = 9 * 0). The corner 


points of this element are 

f>„ » ?•*?!*%* P» 

P- * % h 

P- = Pc - P, » P, - 

P- * Pc - P. - Pc C p, 

The inverse relation is 


(f *♦ i , 7.--1) 
ft • - 1 . ) 

( 1 ’- 1 , 9 - 1 ) ,1 ’ 4t> 


p. « 4 ( %. * p.. - p. f p.. ) 

P> * t - ( R* ♦ P- - P. - jj. ) 

P> * ir ( P.. - p„ „ g, - p. ) 

* ~zr ( P.. - B. - ♦ p . ) 

Note that the four boundaries of the element ( ^ ; + / 
are straight lines given by 


U. 45) 

9*: 1 ) 


p » ( Pc - p. ) ’ (K*V) 

p * ( ?, -p, )* ( p.- p.)p 
p * ( r< * p.)* ( p, •» pW? 
p * ( a -Pc), (p, - pJ? 


ORIGINAL PAGE IS 
OF POOR QUALITY 


(1.46) 


Next, assume that the surface of the aircraft is divided 
into curved quadrilateral elements with four corner points p ++ , 

p + _, p_ + , p . Then, as already mentioned, these elements can 

be replaced by the hyperboloidal element (described above) 

determined by the four corner points p ++ , p + _, P-+' P (see 

Fig. 2). It may be noted that the surface is continuous since 
adjacent elements have in common the straight line connecting 
the two common corner points . It may be noted also that the 
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P C is the centroid of the hyperboloidal element Z^ and hence 
it will be indicated as 

(1.47) 


Pc 


- 

p 


1.5 Surface Geometry for Hyperboloidal Elements 

Next note that the geometry of the hyperboloidal element 
is a particular case of the general equation for a surface in 
a three-dimensional Euclidean space, which is given by 

P = P (?' *?) (1.48) 

where and rj are the generalized curvilinear coordinates. 

Then the two base_vectors 0, , ^ are given by (Fig. 3) 

4 * ■ p. 

% • if » P..-JR 

The unit normal to the surface is given by 


(1.49) 


n * 




/ 4, < 4* I 

and is directed according to the right-hand rule (Fig. 3). 
The surface elemental is given by (Fig. 3) 


(1.50) 


(1.51) 


1.6 Expressions for b^ and c^^ for Hyperboloidal Element 


Combining Eqs. (1.10), (1.50) and (1.51) yields 

/ i 

5 " W$ -tst- *1*? 


(1.5z) 


■I 


H 
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Similarly, combining Eqs. (1.8), (1.13) and (1.14) 


yields 



(1.53) 


These expressions are evaluated under the hypothesis 
that the surface element is a portion of a hyperboloid, 
in Sections 2 and 3, respectively. 


1.7 Trapezoidal Planar Element 

In order to facilitate the evaluation of the co- 
efficients, the integrals are first evaluated for 
trapezoidal planar element (the results are then 
verified to be valid for a general hyperboloidal 
element) . The trapezoidal planar element is a particular 
use of the hyperboloidal element and is obtained from 
this last one by assuming that the two edges c ± I 
are parallel.- This implies that 


p * k * 

(1.54) 

pi -- “ 

(1.55) 

where * is the unit vector in the direction of 

p, • i • e • 


(1.56) 


Note that this implies that 




(1.57) 
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Furthermore in order to avoid crossing of the element 
boundaries, one must have 


Ifl * X 


( 1 . 58 ) 


where the equality sign corresponds to triangular element. 
It may be worth noting that 


h*ih t) u 

which implies 

h'i 1 /,* 


( 1 . 59 ) 


( 1 . 60 ) 
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SECTION 2 
SOURCE ELEMENT 


2 . 1 The Case u o u = 1 

Consider the source integral, Eq. (1.38), for a trape- 
zoidal planar element (S ub sectio n 1. 7^ 

^ ■ ism- 


H -I 


\ I 

I 5 * P.l )( II 5.11 7=== df irj 


-) -I 


(2.1) 


For simplicity, it is assumed that vt = \ /ii*, H Such fnat 

& & U - i (2.2) 

The case u©u = -1 is discussed in Subsections 2.4 and 2.5. 
Note that 


i ^ IjTjT ♦ f zi l 

® 3 | u • & | 


= ^ ^ I + %ou\ - ilnl-f. XU* t*u) J 


I / «, - ^ 

- f= ~r _ ■ _ f===~ + ( 2 , © W J -r 0 


(2.3) 


< ^©a, + /IfliKJ f7| I 

/FFF* /FT 


-14- 


Note also that, according to the first superrule, if - q 

(that is , on the Mach cone ) , 


- \ x 


J. x a x o y * a, * - if & a,) 


(2.4) 


Hence , 


\j±t zJjlRI 

- y* u o £ * u 


= Ha 1 ^ 0 


(2.5) 


along the portion of the contour on the Mach cone (curves M^ 
and M 2 in Fig. 4). Thus, no contribution comes from the 
portions M^ hence 


US, II —5 s 

/FT 


(h ** 


•JT 7 I + 


i ilfxCLU 


- 'if- • 

lillj 

n-'i (2 . < 


since H = 0 on the portion of the lines ^ = £ 1. Hence, 


with 


s * is C K u) ~ V-<>) 

i<(?) * I «* M *1 


(2.7) 


(2.8) 


where H = 0 outside the Mach forecone. It may be worth noting 


that 


r - - 


(y° l)(u° <*)-(% ou) = yxuo y*u < O 


(2.9) 


For, q x u is pointed outside the Mach cone, since q is pointed 
inside it. Hence, Jfoj. < j|® uj and thus. 
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Jf® l * %0U. I O f.r jsMjO 


( 2 . 10 ) 


2.2 Indefinite Integration Along 7 
Consider 


1 • f TF# 1 


( 2 . 11 ) 


Integrating by parts yields 


I * (7-7.) 0 - 1 n~°T t(» 

U - £ o I * 0 




(2.12) 


Noting that 


£ * V-h 


(2.13) 


with 


s ( f* 4 ? f J* 7# (?x+ t f 9 ) 


(2.14) 


yields# applying the first super-rule, 


I(il - (7 -7,)i„ Ml l l il*l 

(If * « II 


I 


■ft 


= 'Hr? 


-+i oi i/f5J * / pj«M J 1 

„ (7 i (hi * i oa) + I^if 5 ] d 7 

= f [ S^M. (ilk t ou)-r i i^ k“ l d, 

)f ioi-<i*or \JfTi h 1 nj'i'i J / 


C jo J. -i»i ),»« 


-« j> o r 


<*7 


j»j t.ou - teujojr. d7 


-II j * Mil 


/R 


j Ill'Ml 1 
\±U°hii iL 

jj.iiil ,/jTf 
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(2. 15) 


Next set 


f. * f/ * ?; w 

ji ^ 

* "super-orthogonal" to u and aj» that is. 


(2.16) 
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This yields 

( d,<- 1. d> - f, T u) o u » o 

( a, ♦ 7. h - f J « ) o d, « « 

or 

£ # T Ci © u - 7* CL © u » fl # © w 
u « ^ - % 5* • 4* - 5. o i* 


which is a system of two equations with two unknowns, 
solution 


h T 
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0,® U GjOPj- q 2 fl >5 

U o u 0,© i- w © q a 0 a ©u 
A, o « u © 4* - uou 


5 , 1 $* oft* fi,. 

(A « 4* OU*i 


0 , tu o 6 t tti 
U * Q U * Ci> 


Note that 


7-7, 


( 7 -(■&»)] *Cj©iL»u 

U < 5* © M « 


J. * Ci ® < u 

-|lu*M x 


Combining Eqs. (2.15) and (2.16) yields 


1 


I 



-II f » u #* 


j7 


(2.18) 

(2.19) 

with 

f* d, ou^t 

( 2 . 20 ) 
(2.21) 


. ( \‘ti» f [ ' in 

) i }« «<“ JpT '* '/FT 

/ l»i b?» u-j»i } “4 *9 „r f i 
i I }» u 

- y - g f £ o u 

i ' e>h Jlf^F JpT r n l 7Fl 


/FI iffij 
*9 . ,r r i . 


(2.22) 


since (see Eq. C.4) 


f/e i * -- ?/e j; 


(2.23) 


As shown in Appendix B, by integrating Eq. (2.22), one obtains 
(Eqs. B. 14 , B. 15, B.16 and B.29) 



f, T - - V" 8 it! ) 


(2.24) 


where 

F(?) = -J— U 

^•4 II £ X J* II 

. /J?gf 

% ok 


Q x ®i >o 


0*0 i =■ o 




hi 
i f*m 


Q x o a* < o 


(2.25) 


On the other hand, as shown in Appendix C 



|? -a, > 4,|* 

-||i* ill* 


(2.26) 


Hence, combining Eqs. (2.11), (2.15), (2.20), (2.21), (2.24), 
and (2.25) yields 


i- ^-1.) t JihihL-i 

I J.* an 

* f. TF 'T' 


-hll hi Xa*C' 

!}'*•'**! nil I}-*, ’ll 


l_ 

-tG*aj 



u o i mu) n n L tHH i hiA 
• ii i * M I) 


♦ fc •)) 





ziiiL zAiit \ l 


(2.27) 


where the principal value, 
is defined as 



of the function 


f 


(%) < \ 


(2.27a) 
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2.3 Source Integral 

Combining Eqs. (2.8) and (2.9) yields 

0 - I* j (2.28) 

where, (note that | M *P »\ * ) if H = l(that is if 

the corner point is inside the Mach forecone) 


T I Si* 2*| 

•U ^5/9/ - -T a TJs a 


-II 5, « 2* fl* 


I r I K f. ■<?,« 


+ (f «4® «,«&) F<7) 


p '//FT I) 


op Z 2 Wal page is 

0P Po °R quality 


(2.29) 


with F(y ) given by Eq. (2.25). 

On the other hand, if H » 0, this means that the corner 
point is outside the Mach cone and there are two possibilities: 
the whole segment -I i "J i. J is outside the Mach forecone or a 
portion of the segment is inside it. In the first case, A Hho 
and hence it is legitimate to set 

X s *0 

In the second case, the integration must be restricted to the 
portion of the segment with H » 1. Hence, the value of I g must 



be evaluated at this extreme of integration/ which corresponds 
to the point where (feq « 0. As mentioned above/ for $eq = 0, 
the first logarithmic term is equal to zero. In the second 
term, F ( ) assumes the value 


f(1) = O 


— w* ( i ® *0 

IM 7 * 7 



5* © < o 

(2.30) 


It may be noted that for q ©q = 0 

£ © &*. *■ f ** A, * / J * | /^ / • 

where N Q is the inward unit normal *vfn to the Mach 
forecone (see Fig. 5 ). Hence, noting that S* is in any case 
directed from ^ * -1 to ^ * 1, it turns out that 


p O 3 i >0 


(2.31) 


if a 2 is directed inside the Mach cone, i.e., for lower extreme 
of integration and viceversa; hence, 

Pf* A® J * - S-*) (2.32) 

with Js-I (>J* + i) for the lower (upper) limit of integration. 
Finally, the last term assumes the value 


- | f • A * A | (- ( £ * A ® \ * &>) ) 


(2.33) 


However, according to the first super-rule for q©q « 0, 
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On the other hand 

$®^ >0 (2.35) 

if is directed inside the Mach cone (and viceversa) . This 
is true if the x-component of u, u x , is negative (positive) 

(Fig. 5 ) since & 0 W >o by assumption. Hence, 

~ U 1 (2.36) 

Finally, combining Eqs. (2.30) through (2.36), one obtains, 
for the case in which portion of the segment is inside the 
Mach forecone, but the corner is outside it, 



(2.39) 
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if the segment — | < rj ± | is completely outside the Mach cone. 
On the other hand, if a portion of the segment is inside the 
Mach cone, then two cases are possible: a) the corner ^ s * I 

is inside the Mach forecone (H =1) or, b) the corner fj ^ i \ 
is not inside the Mach forecone (H = 0). In the first case 
q <j> q > 0 and according to Eq. (2.29), 




I £» * Bj 1 

-M' r 


pios^a,— 


i 


. ( l* a.® a,* a.) F (j) 


- f-isk 



\l 

ffii fUl )I 


with (see Eq. 2.25) 


(2.40) 


Ft?) « 


i i \*i> 

' i i \ t aj 


> o 


i?» 

j. o 5 * 5* <* ft * • 


i 

l&l 



<p 
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(2.41) 
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On the other hand, in the second case, which corresponds to 
qoq £0, according to Eq. (2.37a) 


1 sUo) 


i a, » gj 

-A a, * ft* 


j. * o Fj?) 




(2.42) 


where 


F. ( 9 ) = 0 


a x o a x >0 


1 

h>n 



Q x oa x <0 


(2.43) 


2.4 Planar Quadrilateral Element Internal to Mach Cone 

In this Subsection, it is shown how the results obtained 
thus far can be extended to quadrilateral planar elements. In 
order to do this it will be shown that the second mixed de- 
rivative of the function 



(2.44) 
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is proportional to the integrand of the source integral. 
In Eq. (2.44), 



_J i (lfil II g, II |o Q, 

ll 0,1 1 |1 J. < II 


( 3,0 5 , 


jjj 


Q,o d 


( 

B a , ii 



5, o (?, 


F. ff, 7) = 


l&J I 




H y- * a x II 




> o 
. - 0 
< 0 

(2.43) 

l >0 

K - o 
?, < 0 


(2.46) 
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Note that Eqs. (2.44), (2.45) and (2.46) reduce to the formu- 
lation derived in the preceeding Subsection^ if a^ ® a^ 0. 
Equation (2.44) may be revrritten as 




( i * a,r 


}* <?, 0 


(L * 

1 9, *il 


F. H. l) 


. f. Sl lSl- t^' (ztelk± ) ( 

> la , « <?,/ p (/;^« )•<?.««»/ J 



f * 0,0 ft 


F, ft y) * } « 4 © n F* F?. 




since 


n 


g, » & 

I fli * Q* I 


(2.47) 


(2.48) 
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Next consider the second mixed derivative of Eq. (2.47). 
Note that 


1± 


a, 


t 


ll 

'd r ) 


Q, 


(2.49) 


a 4> _ a 

Q -f 7 

and 


(2.50) 


a £i 

* 7 

Note also that. 


» -f * fi 

for a planar element/ 


(2.51) 

the unit normal n is 


independent of ^ and p 

? 0 s n 

Tf Tf * 0 

Hence 

^ j ‘ «. 0 " j * («■•*>•*)* o 

and 

_L a f nj = o 


(2.52) 


(2.53) 


(2.54) 


Furthermore/ as shown in Appendix B (see Eq. B.5, B.10 and B.12) 




Jf EJ JT°I 

I } « l ii 



{ O Q x > o 




$rj 


nu 
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*7 I 


ll a* H 


fi 


r\ 


t k°Ji 


(Qj © 4 < 0 j 


/if® l 


or 


I 


a') II p II 


(•Per a e O Q, £ 0 J 


Similarly, interchanging indices 

*F. 1 


S'? H l « 


(W a, a a, \ o j 


Hence 


B 2- ( f« a, on , 1 

•? 1 ' JFT 


( 2 . 55 ) 


( 2 . 56 ) 


( 2 . 57 ) 


l *?>°n 1 , — }°^ 

+ . ■ - . J' 4 ,© * 


«)« 


up 11 




( 2 . 58 ) 
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and similarly, interchanging indices 


2 




a, *Q, on 




J-0 Oj 

(f° l)' A 


Finally, as shown in Appendix D, 



(2.59) 



If I* 


(2.60) 


Finally, combining Eos. (2.47), (2.58), (2.59) and (2.60) 
yields (note that fnil *s - fio n ) 


£ls = 1 

non* 



l 


- } 


o n 


tlfl 


- J,x4j.on 
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-5.5 1 'All 

* II J ll* 


i 


C 


Ron 


2 / <?, * Q x ) non - i - 


n l *• * To 


4,o/) }o Qi - yQ^on j.o Q t - j.-n j.Q,i ^ 


(2.61) 


On the other hand, according to Eq. (A. 16), 


- £/ (pti 0 " j 0 *>- P 4>0n J.0 4, - J.-* 2* 5, « *>) 

X 

5 ^ 4 | © Q x x Q t ^.0 Qi %.* Ox o Q, t Qj. J.O 4, + ( J • a, * 


= J.© £ a t t Q>o Q,t 4, 

= f 0 J non I Q, » fix I 


(2.62) 


Hence 


2 lii- . — I fiho h ILiLL. 

Vi**} non \ gl a 


Up II 


_ - \l 

n o n — _ 

//;// 


«• 

' e ii v 


63) 


or 


^ I Q> * Q>] 
9 «p 


(2.64) 
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Hence/ according to Eq. (2.1), the source integral is given 
by (note that for elements internal to the Mach cone, H = 1) 



fTT5 = 1 ,( 1 . 


with 1, given by Eq. (2.44). 


2.5 General Quadrilateral Planar Element 

In order to extend the results to a planar quadrilateral 
element intersected by the Mach forecone, it is convenient to 
use the theory of distribution by Schwartz (Reference 8). 

, For, note that, according to Eqs. (2.61) and (2.64) 



(2.67) 

But while the left-hand side has an integrable singularity, each 
one of the three terms in brackets on the right-hand side is 
nonintegrable within the framework of the theory of real func- 
tion. However, if the right-hand side is treated as generalized 
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function (or distribution) , then the integral of the right- 
hand side (in the sense of the theory of distributions) equals 
the desired integral of the left-hand side (in the sense of the 
theory of real functions). 

It is obvious that if the element is fully within the Mach 
foreconc, then the results described here reduce to the one 
given in the preceding subsection. Also, if the element is 
completely outside the Mach forecone, the value of the integral 
is zero. Therefore, the results presented here complete the 
formulation for a general quadrilateral planar element. 

In order to simplify the derivation of the results, consider 
each of the three terms in the brackets in Eq. (2.67) independently. 
This yields, according to Eq. (2.1) 

I I 

2ffS = (\ J± I 5,- \\ i-jdn = (5, . S.‘ S,) (2.68) 

J, *V 

where, according to Eq. (2.67) 


5 


T 

I 



Qi * 3, o t) 

/!?<) 


If * J y (2. 69) 


5 . ■ - #({ £ ( -Tir • j-itF- “• 


70) 


-I 








(2.71) 


-i ~i 
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where Pf j indicates the finite part (partie finie de Hadamard) , 
of the integral in agreement with the theory of distributions 
(Ref. 7, p. 38-43). Note in particular that (Ref. 7, Eq. II, 

2; 26 with m = - 3/2 and ^ s 0 for x > a) 

Pf f X * (iff f (V * < f(°)— y - '\ (2.72) 

I, Sg ~ ~2 

Hence, one obtains for s Q - L fl 










1 

V 

\k J 


/ 



ffo #>- 

fa ft ) 


ICa) 

fa~ u-73) 

In other words, the singular contribution disappears and should 
not be taken into account. It may be worth noting that the use 
of the distribution theory, in particular the finite part of 
the integrals, is fully legitimate, because the sum of the three 
integrals considered in the following is a regular integral. 

Hence, the three singular contributions (which are not taken 
into account in the theory of distribution) would cancel each 
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other in the theory of regular function, if the integral is 
evaluated in the limit (starting from a region of integration 
without singularity in it and letting the region approach the 
one under consideration here) . The use of theory of distribution 
however is preferable since it yields the same results with much 
less complication. 

Hence, in order to evaluate S, consider first S^. Integrating 
with respect to over the portion of the element with H * 1 
(Mach precone) one obtains (see Eq. 2.58) 

I I 



(2.74) 

where, according to Eq. (2.73), 

H - 0 outside the Mach forecone 

& \ inside the Mach forecone (2.75) 

In other words, the portions along the intersection of the element 
with the Mach forecone (lines in Fig. 3) yield no contribu- 
tions to the integral in agreement with Eq. (2.73). It may be 
worth noting that 

JJ f if m K ( j- %) 


(2.76) 
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where ^ is the value of for which I £ d « 0, and )C is a 
suitable constant. This justifies the use of Eq. (2.73) in 
Eq. (2.74). 

In order to discuss the integration along , it is con- 
venient to rewrite Eq. (2.74) as 

5 , - 5 * - 5 -.' (2 - 

with 



Next, consider the two following possibilities. First 


(2.78) 

(2.79) 


5 0 < 2 ’ 80 > 

1 


(2.81) 


is completely outside the Mach forecone. In this case, 

s: * o 

* 4 

In the second case, introduce the abscissas ^and such 

that H. I ) . 0 -f - 7,* 

-- 1 ?,•<?<? 
s 0 - 1 (2.82) 


In particular, * -I , = i if the corner I ; fcl, y-ij 

is inside the Mach forecone. Using Eq. (2.82), Eq. (2.78) may 
* be written as (see also Eq. 2.53). 



m 


^ - “liiir 

A: ' 

or, according to Eg. (2.58) 


(2.83) 


where 


9, * - ^ ft a, on ( F, l)- F, kf* > ?)J J 

= x,, ft, i; - is, (-'■ 1 > 

* "sTa (M o;r J 7s , F 'fc*' ’) 
I y,f’ 1 ' 1 ) ' “ ( f ' 5 ' 0 ” V ' Fl 


(2.84) 


(2.85) 


( 2 . 86 ) 


Note that if 1 (i.e., if the corner (l,l) is inside the Mach 

forecone) , then-£^is the value of at which H, becomes equal 
to zero, i.e., the edge ( crosses the surface of the Mach 


forecone . Hence 


iif i ‘ 


^ % - %' - 7- 


(2.87) 


This .implies, according to the first superrule, | Jotf, j = 11 J* 
and therefore, according to Eq. (2.45) 


4 X QQ> >0 


a,o 0,^0 (2.88) 

*Here, it is assumed that ^ » a^ ^ 0, that is the centroid of 
any element is not contained in the line =1 (— )• 
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Similarly 

f, (V, ') s 0 a,°«, 20 

_ _ <2-89) 

* - ia7 a ^ ° V 

if 1 . in order to obtain the expression for sign 

Cqoa^) the same reasoning used in deriving Eq. (2.32) may be 
employed. For, note that, if q © q = 0, the vector q is tangent 

to the Mach forecone (Figs. 1 and 5). Therefore q c (see Eq. A. 5) 
is directed like the inward normal to the surface of the Mach 
forecone at P. Hence 

f ° 4 - | C . O, > O (2.90) 

is inwardly directed and vice versa. Note that the case q © a t = 0 
implies that "a^ is not tangent to the Mach forecone (see Fig. 1) . 
On the other. hand, a^ is necessarily directed from = -1 to 
= +1 since along J = 1 

4 * f<- % * i ?,.) (2.91) 

Therefore, if a - i > a^ is inwardly directed at "f* (lower 
limit of integration, see Fig. 6), while, if 'f. f 4 ”1 / a x is 
outwardly directed at ^ (upper limit of integration) . In other 
words 

ft/* (fo5,) * 4 1 

■=--) (2 92) 
Similar results may be obtained for 5^. In summary, it is 
possible to rewrite the results in a more compact form, as follows 


s, = i 5 , (i,i)- ir, a-')- h, r- uj- 1 


(2.93) 
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where 


of 


= O 


i u (t <)- is s 

i 3 «.-«) - i,,a-')so 


(2.94) 

(2.95) 


if the edge /) = 1, or fj = -1, respectively, is completely 
outside the Mach forecone. Otherwise 

V 1 -. ^-*73 


96) 


if the corner = S^, ^ = S 2 is inside the Mach forecone, while 

(see Egs. 2.88, 2.89, and 2.92) 


^~5 t 5+) — 




5,0 Si 


< o 


/JOfl 'As 

if the corner is outside the Mach forecone, i.e., if 




< o 


(2.97) 


(2.98) 


(.2.99) 


In writing Eq. (2.97), Eq. (2.92) has been modified to read 

5*i y* (J.oQ t ) * + i S’, = - 1 

C — 'f 'fpf S, = -f- I 

Next, consider the second contribution, S 2 , to the source 
integral. This is given by Eq. (2.70). Interchanging the order 
of integration and repeating the procedure used for S^, one 



, — L p,("( haj>r\ in 

non ' 1, ^ " nt* /'fs-t 1 


(2.100) 
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where 



outside the Mach forecone 
inside the Mach forecone 

( 2 . 101 ) 


or 


H. (t I, yl- o 

= 1 


-7 i 7, 1 
7, 1 < 7 < 7 


4 

X 


- o lx ' - 7 - ^ 

Hence, Eq. (2.100) may be rewritten as 

^ . 

where 

A ' 

on) -4- ^ 

non k.’ V if* 


s / - sv 


S. ’ a. a f 111 




Note that 


5.- s o 


( 2 . 102 ) 


(2.103) 


(2.104) 


(2.105) 


if Hq = 0 along the edge ^=11, i.e., if the edge fj- + 1 is 
completely outside the Mach forecone. Otherwise, (see Eq. 2.59) 


' 1 S , < 2 - 106 ) 

where 

0 7 /-) 

If. (l I,' 1 ) ~ 4 * n ]^.*i f»(* *' 7/'J 


(2.107) 
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Note that if ^ 1 (or J) - * -1) , 

Xy.fS,,^)=0 

= J-= i-aon- 4 - 4 Si/niioSA ^0«<o(2.108) 

jfon f aa 2 u 2 tr / 

where 

s] r (f o5j ) = -^ = *■ 1 5-*=- 1 

- - 1 5“* * ■* 1 (2.109) 

In summary, it is possible to rewrite the results in a more compact 
form as follows 


* ly, (I. 0-^(1, -0- **.(-1'-') (2- 


110 ) 


where 


or 


I,. (I, l)~ I r> - o 

lrj-u)- if.r-i.-'J- 0 


( 2 . 111 ) 


if the edge = 1 or *^ = -1, respectively, is completely outside 
the Mach forecone. Otherwise, 


‘ja } <5 '°' n U F ‘ fy - rJ ( 


(2.112) 


if the corner = S^, S 2 is inside the Mach forecone, while 


l sJs,.sA = 0 


- 1 

710 h 


O x 0^20 

f i*Q t oh ] — rr, 5^ Q,OG,^ O 

If-S'IM a * (2.113) 


if the corner is outside the Mach forecone . 

Next consider the last contribution, S^, to the source inte- 
gral. This is given by Eq. (2.70). Note that for planar elements 

j_ • h = 


Const. 


(2.114) 
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Therefore 


s 


3 


111 

nor) 





= -1 1 1 2. VT> (2.115) 

non 

with D as the doublet integral given by Eq. (3.28). The evaluation 
of D is discussed in Subsection 3.4 under the less restrictive 
hypothesis of hyperboloidal (nonnecessarily planar) element. The 
results obtained there are applicable here and this, according to 
Eq. (3.39) 


n on k 


(2.116) 


where I D (S^, S 2 ) is given by Eqs. (3.40) to (3.43). 

In conclusion by combining Eqs. (2.68), (2.93), (2.110), 
and (2.116), one obtains 

2TT5 - l s IsH* >)" IsC-t'-O ( 2 i: 


with 


' n on 


where .1^ (5,,$*Jis given by Eqs. (2.94), to (2.97), (shSijis 

given by Eqs. (2.111) to (2.113) and Xp is given by Eqs. 


(3.40) to (3.43). 



42- 


SECTION 3 
DOUBLET ELEMENT 


3 . 1 Introduction 

Consider the doublet integral , Eq. (1.39) 

p = < 3 - 3 > 

As mentioned in Section 1, in order to avoid evaluation of 
finite parts of integrals, it is convenient to use Eq. (1.40) 

P--3Sr (3 - 2 > 


Hence, according to Eq. (2.32) 

D -+3i (*.«) - 

with 

K,(f) = 4^0 V. i s H) 

or, according to Eq. (2.38) 
where 

Ip (-5. 7) * +('?«'».) I. Of,?) 

with I given by Eqs. (2.40) to (2.43). 

& 


(3.3) 


(3.4) 


(3.5) 


(3.6) 
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3.2 Evaluation of I Q 

Consider first Eq. (2.40). In order to evaluate I Q , it 
is convenient to use the definition of the unit normal N 


N 


^ 

I ci, * I 


(3.7) 


and recast Eq. (2.40) as 


I, Cf. i) - 


’ . f-f « a,e ij -J - In 
£)ori l r li, 1 ' 


r. I 


I if* 5,1 


+ F (j) 

- l-H W W 

p \ Ij# / y 

(3.8) 

Note that Pq appears only in q = P - Pq, (Eq. 1.27). Hence 

(N © * o 

(w 8 'ji 


(/J 8 *,) M 


*i7j' 

where N c is the conormal, Eq. (1.15). 
to a^ and (Eq. 3.7). Thus 

sl c o a, i • ci , s. o 

N‘# 5 Sj • z. O 



Note that N is orthogonal 


(3.10) 
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Thus 


K) 0 V, ( <?aJ * - a/S tf, O ~ x 

u l a>* &*\ 

- - ( N°© a, d,o a x - a/^ <3. © d, a,* = 0 

Mo7, (itd^oH) * - w s * - vs 5* g A * 

s - ( jj c ® a, Q x oa x - dt,6 Qt. d» • ) /\ & * &l 50 


and 


MO V. f J. ‘ A// - - A/ c - ^ - A/ © V 


(3.11) 


(3.12) 


On the other hand, for a ^ O a^ > 0, 

/Tf®T /IZ®3, f ©A, 


s * '• (dr *•(“ 


)) 


8 f* III 

= 7 xr * eV -{ 




1 

r » 

f-fai 

Iff,*] 

'Jj + 

im 1 

[ 


(-a/© \) j 


\ 


IA,l 


a, •*,-(}•*)* J J7f 


A,oj,+ 


bl c o$. a.oa, - V c <?<?, Jo 4 , 


-I £* *, t 



| a.ofl, 

-II 


d'l /Ja.otf, &} tf,<>4 *) 

-nj^.ir/P7 + - 1 f * 4 ,h' ) 



f • ji !_ 

-II J* 3,1* J|7j 


Furthermore, by setting 

5* s Si]* ( I ' ^ ) 
one obtains 

N ® V. 4^-) 

III?- II 1$-^.*^)/ 


(3.13) 


(3.14) 


S 

n ^ t (m>« g>r 

\o\(\-a^hy 


% 


* ('!)* 


- V e « fl, ® j.* 3, q(- ti c ) * fl x 

ffo \ S* (\ • o, * 




>A/ C < 9 ? 1 

s*('y*A) 



— i — : — Tu 


x) 


fei(pA,‘Zj+ (yio \*&xY Kp>(i--*y*) 


Tsi X 
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( n fc * a, ® f / 5* ♦ <?, « ^ s a A ) £ ® j 


f • 


- J.*4© £ ya t *a x + l*jN'- 4,* a x 

flfcttiWAL PAGE IS 
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In order to obtain a simpler expression# note that according 
to the first and second super-rules and Eq. (3.10)# 



(Fl c *a,o \*&2 ♦ J.® $• (fr-a,*a 2 ) 

- (J c 0f 5,©^ ~Aj c oa x a,o'p r \o 7l c a,oa x - \oa x a,OA/J) 
x ( \ o \) ( \ * 5, x 5*) 

- (£<?v 0 \-a l *a x + \o\ y c ‘4,* a x )(i.ol «,oa x ) 

+ [\04* \ 9 °>) 

= (}‘Nf 5,#ij) f Cf •«,*&). (f®£) (a,o*>)(M t ' t 4* *&j) 

•t G &i J.® Q+ A/ 4t *4* t 4, ^f0 4 X ^ O *J • Q ,*.4 t 

— \'*i f • 3, * & ( Jo J 4, • fl* -r £ ® *. i-oa*) 

- J. « } sio a,* a x i\o\ a,<s> a x „ ^ 04 , £ 0 


I 


5Y( f# f 5>i * ?* l i°h) 


1 1 * a>\ 


i 

( 0 , * 


J.o \ (a, > a x ©4#* a A )(f-o \ a,oa d - f.o4, foaj 
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m \nul* 0f ioa,) 

• \o a*(4 t ot, \q4 x -\o 4, 4,0 «>)[\o\ \o«, boa>) 

- \o 4,(^0 4, lot, - j.0 4j. & 

- j,o\.(a,*a x o4,*4,){fo^d,oa > ^fs4 x f.o 4>)J 


l Ti kxifo 'tt* 1 oa, a,Q4>-(f,oaj Jj.o4, \oa t 


-(l*5+)*4,Q«, lo\4,&a> -Cfo 4^4,04, j 0 4,- j.04, J.O 4 X ^Q f. (4,0 4 X ) 
+ (lo*») (J 04 J 4 , o a, -(£ 0 */J 4*®^ £© 3 ^, 04 , -(jo fo a* 

+ \oh \ 0 Q X \o\(4,oa x ) - (f.Q4 x f (po 4 , f 4,0 4 X J 

= —■-■ _ — f f*4,(-(fo% a, 04 , - pofl, 4,oa t )( <j.Qf.a x oa x -(j.o4j ] 

I 4, * 4* I / 

( ( £0 4, *l0&>-\04> 4 ,o*>)( f®£ 4,9 4, -($04,)] 

— JjJJ J-f}04Up4,Oa x r4,Kj*a,of*dJ- \Q4jylo4,* 4j}'Zoj*4,)j 


(3.16) 


Furthermore, as shown in Appendix C (Eq. C.14) 


| o%ll-4,ia$+(U*i°\*A>V % i 


T • ■ S — ■ * 


Hence, combining Eqs. (3.15), (3.16) and (3.17), yields 


Mo 7.W ( iIsIl! LA ) 
\ If 1 1 f • «.»$»/ / 


(3.17) 
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2 %«4,oy jo/, J ?«<?»)*- f > 4,oy/o4,JE< 4,1 | 

if**,? if*a t i" 7 J 1 


•f 


(i.i.a -Hi — ili ) 


(3.18) 


Next/ it is shown that for any value of a 2 o a 2 


No V.F(y) * ^*1 




**t f\9 1 


(3.19) 


where F is given by Eq. (2.41). For, if a 2 Q a 2 > 0, repeating 
the same derivation as shown in Eq. (3.13)/ 

Kjor.F- K ln (JHMKdHj ) 




f o a> 


-If** i»‘ fJZJ 


(3.20) 

Similarly, if a 2 © a 2 = 0, then - (foSj * j.xa t o | X S t * -H f*^u‘ 

^ No V.F * NO V. (JHL ) 

\ / 


%• V - 


f **x 


jt°t 


(3.21) 
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Finally, if a 2 © a 2 < 0 






f Q S A 
1 f * a* 1 


-) 



I | \ 

J-K. JoJ if*3j.r 

_5.v_Lfl L 

r Ipial* fj! 


£ 


• fi x o a x 


(3.22) 

Finally, combining Eqs. (3.8), (3.11), (3.12), (3.13), (3.18), 
and (3.19) yields 


J »- * 0 * J * - rsrf' l ' a) * 9 z ( «h *« ) 

+ (t.a.«Ai) Nev.F- S» ( y 7 t*« 

- S'. ( ijeV<y>») jl°L 

\l}llyZ‘Hlj 
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I 

No# 




*(p*) 


\ • I 

Jfof 




tilt [ 

Ip «J' JJoj 


- S„ 1* N r- $ ". f % * 4,0 a Jj-iL- z -Pa J o*l J1 h. 1 

t JJTJ K r r Hftajr) 

- s. c-a» a tiT' ( - }'*•.* !•?•.) I 

v Ilf n j£‘*,**J / J 

* ^ f*» f f tf.) W f - P 4 ' ° f * 4* 

' J7°3 lpj,*2ji 

w;fci±ii* . 


(3.23) 


3.3 Direct Integration (for elements completely inside the 
Mach cone) 

3,3 Quadrilateral Element Internal to Mach Cone 

Consider Eq. (3.1) for elements internal to the Mach forecone, 
In this case/ H S 1 and Eq. (3.1) reduces to 

D= ifj*(iT ji)°* 

*1 "I 

= M{(~^W^ (l ° l)c " l3,rSjl ^ 
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i I 


-I 

2TT 


II 111* 


-I 


f 





since V (fo \) z * ^ . Consider the function 

2, (-(, i) = t^ : li' A* h j 


(3.24) 


(3.25) 


As shown in Appendix D (see Eq. D.17) 

e>^ v*l i J" il f 

Therefore, combining Eqs. (3.24) and (3.26), one obtains 


-2.17 V =■ 



(3.26) 


= 1,(1, l)- 0- h(l.-l)*l,(-L-l) 


(3.27) 


3.4 General Quadrilateral Element 

In order to extend the above results to elements partially 
inside the Mach forecone, it is necessary to use the methods of 
theory of distributions, which were introduced in Section 2. For 

pi 

the conormal derivative of - - is not integrable in the theory 

ip 

of regular functions. The interpretation of this fact is given 


/ 
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in Ref. 1, Appendix H # where the doublet integral is obtained 
in the limit (by replacing _tL with a suitable continuous func- 
tion and then taking the limit as the continuous function approaches 
the original function) . 

In the case of a planar element, the problem can be circum- 
vented in the way used in Subsection 3.2, that is, by replacing 
the integral of the r 'rmal derivative with the normal derivative 
of the integral. On the other hand, in the case of general 
hyperboloidal quadrilateral element, it is still possible to 
evaluate the integral, however only by using the methods of the 
theory of distributions. Consider therefore, Eq. (3.1) which 
for an element inside the Mach forecone, may be written as 
(see Eq. 3.24) 



(3.28) 

Using Eq. (D.17), one obtains 


-2UV - 


( < H - __ * s ' 

) TpHtfp ^ 0 1 <*7 (3 - 


29) 
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where Hq, defined by Eq. (2.75), is used in order to take into 
account the fact that the contributions of the intersection of 
the element with the surface of the Mach forecone is equal to 
zero (in agreement with Eq. 2.73). This is the same procedure 
used in Eq. (2.74). 

Next, introduce the abscissas t » and ^ 

such that 

H* (i 7 } * o -/ - 7 £ 7/ 

- ' 7 " 1 - 

“ 0 rjl ± y £ J (3.30) 

Note that ^ = + I if the point (1,1) is inside the Mach forecone. 
Similar possibilities hold for y* t tj' . Using Eq. (3.30), 

Eq. (3.29) may be rewritten as 

!)~ 

J/j 

^7, 1 1 *4,1 r / 


(3.31) 


or, by using Eq. (D.15) 

’ft! 

— 2TT V 




®I| 

V ?7 

1 ,( 1 , ( 3 . 32 ) 


It should be emphasized that if (i /, “jJ s. o f then z and thus 

1,(11' 9,*) - t (*!,%*) S 0 < 3 - 33 > 


< 
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Note that, if jf. o , then tj. II .» C> and hence, Jt&= s 'j.oq J.o 4 X 
and therefore 


I, ( 1 , 7 .*;= s-r/n 


} 

= <\f»(j.ov,) T'ip'(ioV‘>i/"(j-a,>3 i ) {3>34) 

with sign ( ) given by Eq. ( 2 . 109 ). in order to find 

the values of the other terms, note that, for *^= 1 , and 
are the roots of 

f®f •rf.Tt,)-'9(h* h) 

- (f.*b)*(frh)^fob) tC T*fyj(kh) t( W>)‘ t > (3 _ 35) 

and thus , at ^ = 1 and >7 2 . y * 

Jos, 0 cp, * ) 0.36) 


while 


t'Z’Z J Uf'*f')+ 7/ ft* ?,)]■&’ it f, (, ) <3-37) 


Finally, it may be noted that j.0d,= f f '4i —0 (or, = 

^ - 0 ) at rjs y* implies that a^ (i^ respectively) is tangent 
to the surface of the Mach forecone. Hence, one can conclude that 
J0(7j.= 0 and %.0 4i = 0 cannot occur simultaneously, otherwise 
the element is tangent to the Mach forecone, contrary to the 
hypothesis of small perturbation. Therefore, if l°«! - 0 . 
the configuration is as shown in Fig. 6 and it is evident that 
a more convenient way to evaluate D is by interchanging the 
order of integration, which yields 


-iirp * i. 0- 1. f 1;, 


(3.38) 
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Obviously Eq. (3.38) is equivalent to Eq. (3.32). 

In summary, it is possible to rewrite the results in a 
more compact form, as follows, 

-iTTJ> = l.-O- l >( Zvf-I-') 


where 


ip ( I , l) ~~ Xp (~ I y + 1 ) =- 1 


(3.39) 


(3.40) 


T, - i» (-1,-0 2 0 


(3.41) 


if the edge 1 or ^ = -1, respectively, is completely outside 


the Mach forecone. Otherwise 


I/?, r.) = 

' '• ' up j.i'g 


(3.42) 


if the corner ~ *>, , ^•*•5’*/ is inside the Mach forecone, otherwise 

(*,, r. ) = 0 . 43 ) 

where the term in brackets is evaluated as indicated by Eqs. 

(3.34) to (3.38). 
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SECTION 4 

METHOD OF SOLUTION 


4.1 Introduction 

In the preceding Sections, the problem has been 
formulated and the expressions for the source and doublet 
integrals obtained. The method of solution is 
outlined in this Section. This requires a discussion 
of the value of E on the surface 2. and an analysis 
of the role of the diaphragms. These are presented 
in the following Subsections. Then, the general 
method of solution is presented. For simplicity, super- 
sonic trailing edges are considered so that the contri- 
bution of the wake need not be considered. For subsonic 
leading edges the wake may be included by following the 
same procedure used in Ref. 6. 


4.2 Value of E on 1 

Consider Eqs. (1.6) and (1.8). In order to evaluate 

^ * 

the value of E on the surface, , it is convenient to 
obtain the limit value, as the control point, P„ , approaches 
the surface, Tj , of the body. Following the same proce- 
dure used in Appendix C of Ref. 1, consider a small 
neighborhood Z £ of the point P^ on the body (Fig. 7). 
Equation (1.6) may be rewritten as 


2 ii t 


. <r> 

O i , 



M 

Ilf if 




2 A o c 


\ 


(4.1) 
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where 


ZnSs.&flLl L -4>±(JL)\dZ 
£ y \-dti c Will 3A/*M Ipl'J 


(4.2) 


Assume that the surface 21 £ is a small quadrilateral element, 
as shown in Fig. 8. If £ is sufficiently small, the 
surface may be replaced by a quadrilateral planar element. 
Then, if P 0 approaches the point , the value of ? • l t 
tends to zero; and, neglecting higher order terms in C 
(wliich includes, in particular, the first integral in 
Eq. 4.2) in analogy with the results of Appendix C of 
Ref. 1 / one obtains (see Eg. 3.29) 


2 n Xx. ro. (£j> — — — f — — — ^ 

■0.0 * -O T> f 2 ^ /j ' 


p.-p* 


n-p* z t 


V 


r Aa*** 

* T-P 


* J 


f 1 ' r 

Ho 

tied 
l r 

f' 5 - 

- 

• 

•sf 


v “• 

■iif n 


/./ 


j 





-4, 


Xv W\ 

V4 L 


f 


. a X6I O 4 . 

—i ' : r_ 


R}tl *«L 


‘ l ■-« 


|7-V 
r* *}.’ 


= i » 


(4.3) 


where the upper sign hold if P is outside the surface 
(region E = 1) and vice versa. Hence, for infinitessimal 
values of i , combining Eqs. (4.1) and (4.3), yields, for 
P t °n 2 


2n K h 



. 4> £. _L 

* . W* I! 


•At* 


) 


(4.4) 
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where 

- -o + '■ - > - i ~ 2 ^ ouhxdt t) 

o *4. *4 (?. .>s.vie Z) (4 ‘ 5) 

-L 

in both cares ,'P 4 inside or outside Z ) , in correspondence 
with Eq. (C.13) of Ref. 1. Note that Eq. (4.4) may be 
included as a particular case of Eq. (1.17) by extending 
tie definition of E as 

E = 1 Outside Z. 

=1/2 On Z 

= 0 Inside ~ . 

(4.6) 

Finally, it may be noted that Eq. (4.4) is an inte- 
gral expression relating the value of 3 at one point P 
of the surface Z to the values of \p and o<p/0.\) c on 
the surface Z . The values of / are unknown, while the 
values of .'>•/ 0 are known from the boundary condition, 
Eq. (1.21). (except for terms of the same order of magni- 
tude as the ones neglected in the process of linearizing 
the differential equation for the supersonic potential 
flow) . 


4.3 Numerical Procedure 

As mentioned above only wings with a supersonic 
trailing edge are considered here. In this case, the 
wake does not affect the wing and, thus, it can be 


/ 



ignored.* Consider first a wing with subsonic leading 
edge. In this case Eg. (4.4) is the desired integral 
equation and can be solved as follows: divide 2L into 

small quadrilateral elements, assume 'P and <2 r J /dN c 
constant within the elements. Then, Eq. (4.4) (writ- 
ten at the centroid, P h , of the element Y) yields 



where 

* . \ 9S U.8) 

*• j p p, 

is given by the boundary conditions Eq. (1.21), while, 
approximating the element with a quadrilateral hyper- 
boloidal element, c^k and b hk are given by 


c V -it] 

i 

j Y J 

j r r>] 

(4.9) 

n< [ -Jl J J Z K v 1 Ipl 1 




i 'ijrfl JL ill,, 
jut 

~k k 

] - ’ 

[$] 

*• J j. 5 * p - p 

(4.10) 


with D and S given by Eqs. (3.38) and (2.117) respectively. 

Next, consider a wing with supersonic leading edge. 

In this case, Eq. (4.7) should not be used since the system 

may have a determinant equal to zero. In order to show 
this, consider the case shown in Fig. 9. For simplicity 
assume that the elements are such that the Mach forecones, 


* The contribution of the wake may be included in the 
same way used in Ref. 6 for subsonic flows. 
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C x # and C„ , with vertices in P, and P 2 respectively# intersect 
only the elements 7, # and £ , . Assuming that the element 
*2. is planar # then 


C = O 




(4.11) 


since 1 x # . , £ O on 2- , • Furthermore note that 

* f * 

according to Eq. (F.7) (or Eq. (F.13) with E * 1/2)# 


«=r 

^ c 

h ^ 


« - < s> L ( JL ) c\t - x 6 - 1 (4.i2) 

‘ * j { UjS / » £ 


On the other hand# 

Cm ? * 
r O 


Ks. 2. 

A»-i ... 


(4*13) 

since all the elements except X ^ and are outside 
the Mach forecone C^. Equations (4.11# 4.12 and 4.13) are 
equivalent to 

k z 2 


C ,^ - 1 
c O 


2 


(4.14) 


Similarly 


C i'«~ " • 

r O 


h = / 
H j I 


(4.15) 


Therefore the first two equations in Eq. (4.7) are 

• • 1 < ?, ? . 

I i , r ~ 

. 1 1 j c 9 *) ^ 




t, ( 


(4.16) 


/ 


which shows clearly that the determinant of the system 
in Eq. (4.7) is equal to zero. Therefore in order to 
solve the problem a different procedure is used. 

Consider the Green theorem for the function E 
as defined in Fig. 10. This yields a systems of 
equations similar to Eq. 4.7 where the summation how- 
ever is limited to the upper surface of the wing. 
Similarly for points on the lower surface the summation 
is limited to the elements on the lower side of the 
wing . 

For wings with subsonic leading edges the deter- 
minant is close to zero if the edges are near sonic. In 
this case as well as for wings with leading edge; par- 
tially subsonic and partially supersonic, the solution 
may be obtained by using a diaphragm to separate upper 
and lower sides of the aircraft. For the elements on 
the diaphragm, both o and d'? / j iNl f are unknown, 
while two different equations are obtained by writing 
Eq. (4.7) for the upper and lower side, respectively. 

The solution of the problem is obtained from the system 
derived by writing Eq. (4.7) for the upper and lower 
sides of the body and the diaphragm. 

The method described in this section was used 
to obtain the numerical results presented in Ref. 9. 
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SECTION 5 

CONCLUDING REMARKS 

A general method for solving steady supersonic 
flows around complex aircraft configurations has been 
presented. The extension L> oscillatory flows is 
presented in Appendix E. Numerical results for steady 
flows around wing body configurations and for oscilla- 
tory flows around finite thickness wings are presented 
in Ref. 9 / and indicate that the method, besides being 
intrinsically general and flexible is also accurate 
and fast. 
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APPENDIX A 


SUPERALGEBRA 


A. 1 Super-product 

As mentioned in Section 1, in order to simplify the algebraic 
manipulation for the supersonic flow theory, it is convenient to 
introduce a special algebra, called supersonic vector algebra 
or super-algebra . In addition to the rules of the ordinary 
vector algebra, the super-algebra includes a supersonic dot 
product or super-product 

(X © b - A* b» - Aj by - Af bj 


The additive and distributive rules are obviously valid for 
the super-product. Note that aoa is 

a o & z o hr A 4 i J J7~* "a* < A * 2 > 

that is for a pointed, respectively, inside, on, outside the 
Mach cone (Fig. 1). Hence, in addition to the ordinary norm 
of a vector (or dot-norm ) 

U I * J a - a (a. 3) 

it is convenient to introduce the supersonic norm (or super- 
norm ) 

(an -JilTTi (a 4) 

Finally, it is convenient to introduce the concept of covector 

k .5] 
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With these notations, it is immediately verified that 
d @ F - c * 5 s & * t* 


(A. 6) 


It may be worth noting that 

fl • f * c * Q ® £ j t C - & 0[b * C c - 1 J ]j C * c C (A. 6a) 


A. 2 First Super-rule 

Throughout the subsonic finite-element formulation (Ref. 6) 
the following rule is used 

a $)•(£* d ) - (&• c)(t* <$)-(&• d)(b‘C) (A. 7) 

The corresponding supersonic rule, called for convenience, 
first super-rule , is also valid 

(a ' h)° (t'd) = (a® c)(Eof)-(aod)(toc) 

For 

(a* Do (c* 6) 

~ ( &§ b* - b*) ( C* 6 k • Cm d#) 

“ (Ox by - Ay bx) ( Cx^y -Cy dx) 

— Q.y kg Cy Jf &T, by Cg dy - Ay kg Cg 6/ *■ f>/ Cf 6t 

~ Cg a* -At bg Cg 6 g + &+ b* Cm <Ag * fig bg Cg 6 a 

- Ax by Cm dy - Ay kx C, 6 A t Cm by Cy dx * dy b A C.t dy 


(A. 9) 
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while 


( a ® c)( F o Z)~ (cloZJCS o c) 

( Qt C K - Q.y Cf — Q.q C-t)( fa* dg — by dy - b+ dzj 

- (flxdjt- flpdf bgC* ~ byCf - 

Qig Cm (b/d* - bydy - bfdt) 

- Gy Cy ( hgdx -Vpjff-bzdz) 

- a*C z (b,dM-by4,-y£) 

- a* 6* (yt* - b/c,- t>*ci) 

♦ Gy &) ( b* Cy CeJ 

t (Xz di ( bx Cx ” by Cy -bj^J 


— Qy b 9 C y dz «• by Cidj- dy b* C * dy - t/ C, d« 

** Ct d* - A Ct dt * tjr Cr d* 4 &x btC* (1 a 


- Cg by Cy Sy ~ Qy b I Cy dx * dt by Cy d* ♦ £/ hr Cy dy 


A. 3 Second Super-rule 

A second rule of the super-algebra is 

(a 0 d)(l * c) o (It c)-( 4 ‘l*c) 

_ 5 o c cj® ^ £"* 4) + 

a 0 F (c,*)o (Z*a) 

(A. 11) 

Note that the dot product appears in the triple product. In 

order to prove Eq. (A. 11), consider the regular vector algebra 

rule _ 

a * (* * £)- C (*• t) 

\ A • j 


This yields, for the covector a c , (see Eq. A. 6) 

oS Z) = l (i'-D-cCa'-D 

*l( 5 o c)~ c (aol) 

On the other hand, according to Eqs. (A. 6) and (A. 8) 
A c t (l»c)o 4 e « d* c ) 

ss. a c o a c (£* c ) q CS* c) - ( a ‘ ° d* c-) ] 

= ft © (1 * C.) 0 (l x c) ” ( 4 • l * C~) 

while, according to Eqs. (A. 6) and (A. 8) 


( l ( a qI)- c. (Sol)} o (l (a oc)-c (oiol)] 

= ¥0 Wfl® c) * z lo c a oc aol ~ coc. (aol)* 


(A. 13) 


(A. 14) 
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Combining Eqs. (A. 12), {A. 13) and (A. 14) yields 

ZJoff ({a c)o (£* z) - ( tx c) 

‘ » ( a c * O C a‘* (l* C.JJ 


- ( 1 (aoc)- c (&Ol)]oCb (aoc)- c (aol>)) 


- aoc CBx c. -Q\h * &) + Hot ( c * 5)«(c * a ) 


(A. 15) 


that is the second super-rule, Eq. (A. 11). in particular, 
for a = q, bsa^, c = a^, one obtains 

j.® j l> i© < 3 ,< f 4 1 * 

= 1*1 (l.Q' ol, j) 

V** t 

A. 4 Third Supersonic Rule 

A third useful formula, called the third supersonic rule, 


is 


3 * 4 & b f c a • { * j 

-t- A. * J o b * C 4 *2 ■■£ 

The proof of this rule follows 


# * -po £ * C 3 'j * 2 

O 


(A. 17) 


( 4 * 3 o b* c)(a • f < j)- (3x 1 * c ) (a • j * 2 ) 

= (aol co 3 - a o c b o2 ) a * • $. 

- ( q ol coj - aoc. lo}) d* ? '2 

- fiofc fro 3 a >{■ j - co) a * f * 2 ) 

-aoc (l o 2 a *-f ‘J - I © £ a *{-2) 
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= a • V (c e - 2 a x {. cj - c c -j a * f ‘2 ) 

- a-c'f'b'-d a,f.2) 

= i? ♦ i c [ c% (a* f )]*5 * J - a*c c [ l c x (a tf ))’ 2* j 

* 3- (c'-aHf-j' ))} 

- 1 ' c £ ( (b c . f )(a- 2 *)) - ( i c -aj ( f-hp] 

— ( & • y c c°‘ { — A'C* y c . ~£) & • 3 * 

s ( a o £ c o f - a o c b o -f ) = ( a * f o * c ) Q ■ 2* J 

(A, 18 } 


In particular, for a = b = q, c = f = a^/ d = p^, g = a^, 
the third superrule reduces to 


(}’ h • f* i)[pl - 1) l) (l-k • },) 


(A. 19) 


page is 
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APPENDIX B 
BASIC INTEGRALS 


B.l Integral 1^ 

As mentioned in Section 2, in order to obtain Eq. (2.24) 
from Eq. (2.22) is necessary to evaluate the two integrals 

■■ ’lTpr i; 

>■ • Sih^ak ri > 

The integral I 2 is evaluated in the next subsection. 

In order to obtain the integral in Eq. (B.l), it is 
convenient to <_reat independently the three cases 0 Sj <0* 
Consider first, a 2 © a ^ > 0* Note that 

jrj IsJZ * i si >) 

/SsSs + feZ [ Jfo l * a ' 9Q 

_ J a*Q a x 

Note also that 




l « « » 1 Vf. ♦ * ?, )' (p. » ■ * r,)\ © 1 dr. - rf,) « (f. - ? f, i) 


= 0 


(B. 4' 


Hence 

2—1 (i , Jp° iJiiOti r i°i ) _ Aoa 

27 1 



Note that 




ORIGINAL PAGE IS 
OF POOR QUALITY 

(B. 6 ) 


Hence, it is possible to write 


of S t @a A + | 

I l* 8*1 


(5 x o&t >o) (Bi7) 


Note that 


'x s JU (x +/J X*> ' / 
X = Ht\ ( x + J!c~- 7 ^ 


, X 2 I 


(B, 8 ) 
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Hence 

i\ + a* ‘i 

— £/i /[(fofijj f ~p*djO ■* I fc & J 

I \**A 


- In 


1 1 - 4 II 



and 




Next, consider the case Sj 0 ^ = 0. In this case 


“>) * £i((r-*isf''7 s ') 0S >) 


* £ (if. 'I?.) if,)) - O 


(B. 


Hence 


JJ°7 


1 pea. 


®') loi, JJTf JJZ J 


(B. 


Thus 


i * Jt 0 i 
' \ 


(B. 


Finally, consider the case Sj Q aj < Note that in this case 


* S'; 


-»/ S o 4* 


1 


4>o Q t 




J}***y fT Ti o 5 r i}*sj 
y [ifTW 


Hil> II 


J - p *<L-(po 0*)' 11 ? x & 11 


til® fix 
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In summary 


in ill n A* llfd 

II a* H I \*Ol>.\\ 

= £&(i°A) Co, k' JilM 

n a> H ii l x £ | 

* -Sr (i°*J a x © a x > o 

ii a* a up 5, 0 * 


I, 


J i < 

% o a, 


, Q>Q Q x *o 


i. - 


! Sin ' \lL. t ft o a < 0 

l\a x i\ lipxaj 
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B.2 Integral I 2 


Next consider Eq. (B.2). Setting 

l - Q. + rj a> 

one obtains for u • u = 1 


l 


O u 



\ 9 v =■ j + + 


(B. 14) 

(B. 15) 
(B. 16) 


(B. 17 ) 
(B. 18 ) 


(B. 19 ) 
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j. < U O J * Ll - fo f UO U - (£ o Zt )* ■- < 

- ^ - Q* + i ( a b) p + ( {- bO y 

- C + + €.^* (B. 20) 

with 

A = o (if 
b s a t o a 
~ a, o a, 

P - &, <* 0* 

r = © a x 

C - oL - a* . a.o4 t ~ (q,qu) x * a t * & oO,* 
d r p - ab . a,c>a 2 . a ,0U 5iOU * d,‘UOQj.tiA 

£. s )f - b* • i* (mm) 4 * ofl,' w 

Combining Eqs. (B.2), (B.18), (B.19) and (B.20) (see Eq. A. 1 

of Ref. 6) 



*—r..-p‘-U)J*M->c) _L) 

/jtC- d * | /e-C'd' JJTfJ 


Note that 

l a<z.~bj)y + (ad-bo) 


(B.22) 


= & ( Cj + d) - b (d J + c ) 
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= Q t (X-b')7r{p-ab)] - y ((p-ab) d*j) 

- a + - b ({7* *) 

s 5, o u (&>y+ a,)oQ* - £© u (a»y+ q,)o u 
= u \o a x - & J.o a. 

*-p'Hoa.,& .-//is <B ' 231 

Furthermore, 

e.c-A * £o £ (&>a J )o(C(*a i )+j.oa,( f*uo u*& 2 ) (b. 24) 

+ jf'U (i*4**Q**u) 

For 9 9 ' 

ec - & - ( 3 X * u o a^* &) (fi, *U q 4» * fi) - (d,*3 qa,* u) 

- ( ( tfj( u o uj- (5*© u/J (fj. O d.)(a O Ct)- ( a, o u)J 
- a. o a a>*J a 

* (£. o a.)(&° <U- cl° Z j (i.o mJ- 3,o a 9 (a k o u)+ (q,q u)fa 

' %,)*+ 2 5,0 \ 3,0 (A OjOCi - (Q,oiA) (Q x oIa) X 

while (a* (ijo (utaj *r \ o l % c f « *) o (u * a* 

+ i o m ( \ t i3q l&t * a) — 


(B.25) 
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— {.(fi'+'/A,) o (&.*•) fr)] ou* & ] 

+ 1 (fr + 9fr)0 4] ( (4*? fr)* U o(U* fr)] 

♦ [(&.+? fr) 0 &JC ( fr + l&J * £x°(a*,**)) 

= ( fr o Oj * ^ a, o fr * f) *&*) ( U * fr o m « fr) 

+ ( fro fr «■ 7 fro fr) ( a,* c< o a * fr) 

♦ ( fro fr* *) fr o fr) J ( frt U\ O tA 1 fr) 

«• C dou + 'j fro m) ( Q t A fr O fr , tf) 

♦ ( a t o Ja «. 7 Q x ou) *) (fr* fro & x td) 

— 0., o a, ( C\ 1 1{ fro fr - (do fr) x ) 

♦ fr o fr ( fr o ^ (A o a x - Q, 0 fr M O u ) 

♦ 0 0 o U (Q,Qfr fro M -fl, 0 W 4 • 4 J 

♦ 7 j 4 ° 4 ( ^ © u frQfr-(tAofr) X ] 

+ fr o fr ( fro u M ® fl, - ^ * 0 4 you] 

+ fr O (A { fr o fr fro <A - Q 9 0 A fro frjJ 

= fr o fr fro fr - fro a, (ti O frf+ 1 fro fr a. OU do fr 
- (d t 0 frf - fro a. (a,o u) x 
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(B.26) 
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On the other hand, applying the second super-rule, Eq. (A. 11) 
(with a = q,"E = u, c^a,) 


(i-UtHjfe j.O j. ( U < 4, » (Si it O St 

- f • s (i t tu o Sjtf) 

Thus, comparing Eqs. (B.24) and (B.27) 

cc - l* = u* 

Hence, finally combining Eqs. (B.22), (B.23) and (B.28) 


(B.27) 


(B.28) 


i * .. J- u* 

* I f%°T~ I 


(B. 29) 
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APPENDIX C 
TWO USEFUL FORMULAE 


C.l An expression for r. • h 

In this Appendix, it is proved that 


® fv Wd^dJ * - f J ♦ 4, * *»/ 

(C.l) 

Note that 


f - K • ( 7- 4. 

(C.2) 

and 


M 

» 

* 

4 

4 
^ 1 

(C. 3) 

Hence 


?* ~ W * m lr~ °1' 9*) fc * 

(C. 4) 


or, according to Eqs. ( 2.20 J and (2. it) 


4 4 £ * A* O iA * 5 X - t * o J » (j 

K « > - “ ' — — - — U - zr — 1 ■ a. 

-H U * I -I! tf * 

- ~ j/ll'ij'f * ' 

Id,* M / 

+ Q t * Q*) (c. 5 ) 

Hence 

9* h" ° h" 

- - 4- * * 

= l n d,* dj 11 +(£*<?"« *,* 4 ,) 4 , + 
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+ * J 1| a, * 4JI* ( p 5» o 4 , « aj foa, 

+ 11 i>'*X(y *>•*>’*>) l 9 «* 

+ ( l*4,o£* *J 4t°*xJ (c.6) 

On the other hand, according to the second super-rule 

(f- a,* i) - l o l , a> o i, , i t ) 

- J- o a ; (}' «i ® «, « 

- £ • ( f K *> ® * ®</ 


(C.7) 


Hence 


13 ■ M 2 [f/« f * l\a,<ajl* + (%-a,’5j) 


= f°fll a, * aj' 4,o a, ♦ (f*4,*4>*a,) 

+ ^)| d,*d 2 \\ j"( f*dj 0 4, ydj) t° 4x +( %* 4, O 4y d.) f° *Jj 
+ 2 (f*Q±* 4 ,*dJ(y*> & h* 0 i) 


- - T ^ X 


+ // 0, * ^ //^ f ® £ { *, ”****,• *>) - l°*i (iyK ° *• * **) 


- £ ® ^ «** * jJ 


0lil WNAr „ 
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( f- Q x o a, 

X 

,a x ) 

a,o a, 

.a _ 

* ( j * 4,0 ^ 

til a,- sj‘ 

l* 

Qj_0 4? , * 

Qx) i'Q 4, 

+ (f,*a,oa 2 *a,)jfG 

+ a ( * @ X 

oa,. 

aj( y 

4, O Q x x tfy 

) a,* 4j 

4>® a,. 

a x {\ 

* 4» ® Q, 

* a x a,o a x 

+ p 4,® 

i‘U5,*aj \o a,) 

+ %* 

4, o fix* a, 

( p4,o 4^ a, a x °o. 

+ $r *• A* O 

a,* a> 


+ ii a, a a 2 

tpoa x ) 


* 0 

(C. 8 ) 

since* 

f ° Q t S t .if 5f&' =j*a k oti r Q , c Q t o& l 

•t* £ * ® A* * 42 1 • fl. 

t -t : », Z. • !C •. 2 • ! : (C. 9) 

For 

O <2j C?j * 42# O * fix - £© ® 4?j &K 0 Qt) J 

\ 0 » «L0 ) 

* 

For ordinary algebra, this corresponds to the v/ell known re- 
lation between tensor components *jL *■ Q-.l 


5 , 
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i! 


t i 

Li 


i! 


{ ; 


i ' 


T I 


i! 


IT 


I 



and 

* On © 0\ < On ® $| + y. * Of Q 0* * Oi i 0^ ® (?* 

= J_ © $, Q^o Qt d^o - J. q On On 0 d| fljj • 

■r j® d* 5,0 5, Q| o fl K - Jo Q t f Q |0 fl,)* 

* J- ® 5f® Qj[ - r5,o 5)f^ j 

(C.ll) 

Equation (C.8) is the desired proof of E^.. (C.l), 

C.2 A Second Useful Formula 

In subsonic theory, it was shown that 

^ f '&•*&*) + s l f*Gi l (C. 13) 

Here it is shown that for supersonic flow 

ipo,it (i %.tiX (c.i4) 

For, according to the second super-rule (see also Eq. B.27) 

J • i (f-a,*ci)*r < \*a x o i*a x f 

= % • } ( f • J ( 5 * 4 o *, x <5J- J © £ (a, * a* o 4 , « \) 

- Jo a,(a^a t o 5,* j)) * (Jo j - Jo*« 

= QiOQt-ca.oaJj-fo l faUnfaoi' a x *\ 

d* oa t )~ J© l Jo 5, ('^o^tf.oj-*,© j a % oa^j ^ 
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a,oq x foa, }Q4>+(fo«, 

* Cf® l l 

= (f*a,Q } * 4i) = (%* 4i O f*4t)(f*0* & i'Zj) 

* Ilf II II £*^4* (C« 15) 
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APPENDIX D 

DERIVATIVES OF I Q AND I g3 


D.l Introduction 

In this Appendix, it will be shown that for any hyper - 
boloidal quadrilateral element, the second mixed derivative 
of 


V 


-liL W/l-J-r-f ) (D.l) 

/{-•"I in pi lya,* Qi\< n tp 


is given by 

tTI-P _ jr ' 

~ ii } a* (D * 2) 

while, for any planar quadrilateral element, the second mixed 
derivative of 


•Si 


- ryn i bn' - h. 4 V : ), ?. ; (±hk ± ) 

b / j- <j,< «J/ » ! upyi.e>l 


(D,3) 


is given by 


?’j» . j - f • g. ■ q. 
*?*7 ’ ' j j. t’ 

Note that 


(D.4) 


and 


3 - t St 

~ *>) 


j> <?, _ » 

*7 * 


ft 


J _ 

zt 

*1 m 


4, 

a 


(D.5) 
(D. 6) 

(D. 7) 
(D.8) 
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f« f» p° j- }• 


- f.£* \-fr 


-i 


+ 1 1-' fi • ?' ^ (j-* flj o \* 4 ) J- fa i 


Next note, as shown in Appendix C, Eq. (C.14), 

£ ' £ f \ ' 4 ' £ ) + %*Q*) = 'hTr'QA II f* 4* 


4 



4 ,* ^ 

(O.ll) 


(D.12) 



Moreover, note that (see Eq. A. 19) 


((<?„* «,« }« aj pop- <.p,'t>, o p, l) a,o p] ha,.a> 

"•(( I' ft • «J( l' 4* l)- (}’«,• i'SjCp-Pt'i)! j°j 

= ( ( 4**Jp*ai ® \ \ a,o\) p ® p 

- - x 

- ( J94 |oiJ | ' 4 

- ( (j'4® |»tfj ft) Jo j 

= (~?®j 4«, -(poaj](Jee,)(j.- a,* aj 

+ <$•*** pi ] (p- <w. 1 £°J 

® flf'aXff®*’ ?®| p- 4. > f, ) 


(D.13) 

since 

|o j = - ll Jx a x II 1 (D . 14 ) 

Finally, combining Eqs. (D.ll), (D.12) and (D.13), yields 

zh. » 

**) " ipa,t tip^u'tlpi 
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Next, consider the second mixed derivative, noting that 

^f(p4) = ^r((p.~7i)* ( P.'7A)]= 0 (0.16) 

one obtains 




l y a,t 



r 

= lf7|i* ” t y s ‘'P>~ P a ' P*' A l) 


=s (2 %-a,t A + }o l fa -a, a q yj^s. 


-J-oq, QrQ t *Qx- fo <3, jL» Q,t ft 


||f, 0,1* flp 1 [(J- 0 **) ~ l°f *'® ]■' Q '* a * 



(D.17) 
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D.3 Derivative of I« 

Note that for a quadrilateral planar element q.n is 
constant. For 




(D. 18) 
(D. 19) 


Therefore, using Eqs. (D.l) and (D.17) 






nin 



(D. 20) 


in agreement with Eq . (D . 4 ) . 
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APPENDIX E 


SUPERSONIC OSCILLATORY FLOW 


E.l Integral Equation 

In this Appendix it is shown how the results obtained 
in the main body of this report can be extended to super- 
sonic oscillatory flow. Introducing the variables 


X = Y* £- Z-3. ! ; t il-ul 

hi l J - l (B-l) 

A * 

and the complex potential <+> such that 

f <t> (KXZ)e 

(E.2) 


the integral equation for the subsonic oscillatory flow is 
given by 


i a r 

"♦ = f 2 l3F 



(E.3) 


where Z surrounds 


body and wake. 


C.2 Boundary Condition 

The boundary condition is given by 


V (p 


3*1 u dS 
dt * 


(E.4) 
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or 


<? r? a > ^ b .w* *& * o 

%tz 4 + rt sT 3 ?X JX *X 


(E.5) 


where C and 4> ere such that 


Ux+<f -- OifrX+t) 


(E.6) 


Next assume that the motion of the surface consists of small 
harmonic oscillations around a rest configuration/ that is 


4(W) + 2 X ^ Z K 


JIT 


Then/ setting 


jii 


(E.7) 


(E.8) 


one obtains 



f /SJ 

-v . / 2J1 ^ 

-K,Sol7 

? U * 

*^2 **2 


* if^A , 

^ ' JLT ] 

B ( ±,x 

2X / 


\ <SL7 


•TIT 




. kj 3 + a Oil*' 

ax c>x jx jx dx J 

~ *1? UJIT 1 

, <20 C ! s O 


;jzf 


*X x 


J 


(E.9) 
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Assume that the surface is given in the form 



(Upper surface) 

-fZ-? t Cx.Vj.Z 4 fxAje ,Jir ]=0 

(Lower surface) 

(E.10) 

with 


Zu.l IX ' y )= OUi 

(£•11) 

K ( (*■■) - 0(i ' ] 

(E.12) 

or, in general by Eq. (E.7) with 


$.'±[z-z..t(xy)]*°iv 

(E.13) 

M: - ?£,<. , o (t ) 

dx JX 

(E.14) 

o 

it 

?? 

N 

(E.15) 

and 


<v ci \ 


2. * otA, 

(E.16) 

j£c C(t J ) 

*>x 

(E.17) 


onrcrx'Ar, page is 

ut'ALny 




Assume also 


J 1 . O(i) 


(E.18) 


This implies (see Eqs. E.25 and E.26) that 


' O (s j 


y ^ f ■ v 
<p z O ( c j 


(E. 19) 


(E* 20) 


Neglecting the terms which contain e (which are of order 
4 

t ) and separating the steady from the oscillatory terms, one 
obtains 

- U V r B ,;x X 4 <?X 


_ V & © Rj, 2-7 $5> /' & + J < 25. 


’ A £# w. iiLo 


t £L [ Z22 - + ~ : 


(E.21) 


(E.22) 


Introducing ,'■■ such that 


7" * -<7z/'lX 

<p - (p e 


(E.23) 


Equation (E.22) reduces to 




-V,yA 3U »< * 4 ♦ 1 ^ T 5 - ^ 4 

4 £ ;JZS 4 i2r * At 1 f2d f §Li . i JU1 4> ]< 
i 7\ & *>x w( ^x \ .>/ / 

f i>.l £A ]. o 

*>x <>* r 


. 7 9. 

«✓? 


1 


(E.24) 


105- 


z 

Finally, neglecting terms of order £ in Eq. (E.21) and terms 
3 

of order * in Eg. (E.24), one obtains 


_ t7 $ & V b - - - 
V *H o u t* a 


ih 

3 2X 


* / * \ iJir \X 

V £ g V 6 . _ [ i SS2. 5 * k*' ]c 


(E.25) 


(E.26) 


In particular, for 


SL.il 

M J u “■* O' “ L J 


(E.27) 


(where the upper [lower] sign holds on the upper [lower] surface), 
one obtains 


^ - * -f [*- '-./vi 

(E.28) 


- ; j 

(E.29) 

1 r /A/./ . i A/, 

IV s 1 *1 

1 XV* 3 

(E.30) 


and 

It i gglfj a. 

^AJ C /^-4/ 

where 


d) 


, .» yV /■> A' ? 


.£. -f 


(E.31)* 


(E. 32) 
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C.3 Pressure Coefficient 


The pressure coefficient can be evaluated by using the 
linearized Bernoulli theorem, as 

- _ 2 f 3 + J 

( /v j ~f 3 & X / 

For oscillatory flow, setting 

j k ,slT f i M r -^ x ) 

<?> = tf> e « <p e 


(E. 33) 


(E, 34) 


/■O z 1 ^. / 


(E. 35) 


one obtains 




$,- 2(1 -J 


/ 




. .2 /7j2- 


* \ 

i 
/ 

j 

i. ^ 3 J ‘ £> j 


(l.- d 

\M A J 1 


A 

^ + - 


i 

/ 

, -*** 


A 

t -1 I • t 1 € 

£ L Aj X* J 


/A** 


^ , uixM ^ 4 ,.JV ,,'SlXlX 

2 / e y7 € 

7X 7 J 


B 




- -§ * 1 (f 

- 2_ * £ 

* 2X 


K € 


y 



? - X *73 

t' ' 


(E. 36) 
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APPENDIX F. 

SUPER-SOLID ANGLE 

In this Appendix the concept of supersonic solid angle 
or super-solid-angle, dfL* , is introduced and it is shown that 

s siSl — Q ( P <— j 

T 

-.in (P. ~Z) 

«. 4jJ f P X i 

° (F.l) 

Note that (see, for instance, Ref. 1, Eg. 6.6) 



(F.2) 


where the finite parts of the integrals are understood, dd fl 
is the usual solid angle, while 

JJi' , J V !.fl it SI 

Util i 
t 


(F.3) 


8 - 



is the analogous one for the supersonic flow and is called here 
supersonic solid angle or super-solid-angle. Note that 


(F.4) 


depends upon the direction of j but not upon its magnitude. 
Therefore if P ; is outside the closed surface . then (following 
the same reasoning used for the usual solid angles) 


0 dSl' * C 


(l> otJ’iidL z) (F.5) 


Furthermore note if P 0 approaches P + on 2 . y one obtains. 


according to Eq. (4.3) and (F.5) 


(& dSl' ) s f & dSl' - 2$ i /£ ■ dl) I 

' £*<>1?.+?, ( *Z.t t i L ‘"e 'J (f.i 

- 2M t0 


Therefore 


dSl 1 , 2j} 


(P. ) 


(F.7) 


Similarly, p if P is inside 2 the super-solid -angle is the 
0 f 

one intersected by the Mach forecone on an arbitrary surface, 
for instance the plane X -X* =-a. ( Fm. r. 1 j. 


# w 
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= . (F. 8) 

Equations (F.5, F.6 and F.8) are equivalent to Eq. (F.l). 

It may be noted that according to Eqs. (4. 6) and(F.l) 

£ c j - J- £> dSl 1 (F. 9) 

4/? * 

It should be noted that Eq. (F.3) and (4.6) are valid only if 
the surface, L , has a unique tangent plane on , 

while Eg. (F.9) is valid even without this restriction. The 
proof of this statement is not of interest here, but could 
be obtained by following the procedure used to derive 
Eq. (4.6) by using Eq. (F.5) in order to evaluate 6 L in 
Eq. (4.3). 

It is worth noting that the analogous of Eq. (F.9) for 
subsonic flow is 

£ - I , _L 3 clSl (F.10) 

U ~ 47/ V 

The proof of this equation is similar to Eq. (F.9) and is 
not given here. As for Eq. (F.9), Eq. (F.10) is valid 
even if Z does not have a unique tangent plane in 
Equation (F.9) implies that 

2 ht 

- 5 £ t Jri ‘ - 2 

,1; ‘WAL PAGE IS 
J -‘ POOR QUALITY 


. C 


hK 


l * 


z C - 




fl 


#3 


-J— )«£ 

M/ c i. Ml 


(F. 11) 


jt It i. 


r 


- jA-- 


n 


1 f Z 

i * 




- J? ( 2 /7 a ) ■=. 


that is the sum of the coefficients of Eq. (4.7 ) , where 
E = 1/2, is always equal to two. Note that this result 
is valid even if the surface does not have a unique tangent 
plane in since in this case Eq. (4.7) must be 

modified as 





(F.12) 


Equation F.ll may be used to evaluate E for 
points with slopes discontinuities (such as corners 
of quadrilateral elements or the apex of a cone) as 

£ r - 2 C + t (F.13) 

^ A* 


The above results are valid for subsonic flow as well. 



r 



Pig. F.l. Super-solid - angle 


